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Abstract: A new block spin renormalization group transformation for SU(N) gauge models 
is proposed near the non-trivial fixed point in perturbation theory and thereby the expecta- 
tion values of various Wilson loops on the renormalizcd trajectory near the fixed point are 
explicitly obtained. An improved action is obtained as in a preceding paper and a criterion 
for the scaling behavior of physical quantities is also given. 



1. INTRODUCTION 



In a preceding paper p| (to be referred as [I] ) we have emphasized the significance to improve 
a lattice action for the continuum Hmit and have shown our strategy to do it in renormahzation 
group approach for asymptoticaUy free lattice theories such as 2d nonlinear 0(N) sigma models 
and 4d SU(N) gauge models, most attention has been paid to the former. In particular, the 
correlation function on the renormalized trajectory has been explicitly obtained near the non- 
trivial fixed point by block spin renormahzation group in perturbation theory. An improved 
action has been obtained from the criterion that it is most close to the renormalized trajectory 
by defining a distance from an action to the renormalized trajectory. A criterion has been also 
given for the scaling behavior of physical quantities obtained by Monte Carlo (MC) simulations 
with a given lattice action. 

In this paper we will make a similar analysis for 4d SU(N) gauge models. We will not repeat 
in principle the general discussions given already in [I] to avoid the repetition. Therefore we 
assume that the reader is familiar with [I]. 

The organization of this paper is parallel to that of [I], except for that the sections which 
correspond to sections 8 and 9 of [I] do not exist in this paper: The model is defined in section 
[21 block spin transformation is introduced in section [3l the expectation values of various Wilson 
loops on the renormalized trajectory are obtained in section [Hand the block spin transformation 
is performed for various lattice actions in section [5l An improved action is obtained in section 
[51 implications of our results for MC calculations are discussed in section [7] and the connection 
between the existence of instantons on the lattice and the renormalized trajectory is given in 
section [8l Section [9] is devoted to discussion. In the appendix [A] the explicit form of the 
propagator is given and in the appendix [B] the derivation of eq. (j3.9p is given. 



2. SU(N) LATTICE GAUGE MODEL [2] AND WEAK COUPLING EXPANSION 

We take a four-dimensional hyper-cubic lattice and denote the sites of the lattice by an 
integer- valued four vector n. The link variables of the SU(N) lattice gauge theory are N x N 
unitary matrices and denoted by U^{n) for the oriented link from n to n + ft. For the oriented 
link from n + /t to n we assign U^^{n). 

There are infinitely many choices for lattice actions which give the same naive continuum 
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(c) (d) 
FIG. 1. Four types of Wilson loops in the action. 

limit. For example, we adopt four coupling gauge action, including up to six-link loop: 

S = — {co Tr (simple plaquette loop) + ci Tr (rectangle loop) + C2 Tr (chair- type loop) 

+ C3 Tr (three-dimensional loop) + constant}. (2.1) 

Four type of loop are depicted in FigdJ In sum over loops, each oriented loop appears once. 
In the week coupling expansion the link variable is parameterized as 



U^{n) = e 



iagA^{n) 



(2.2) 



end the action is expanded in terms of the field A^j(n). Requiring that the action (j2.ip reduces 
to 

v2 



(2.3) 



in the naive continuum limit, we obtain the normahzation condition 



Co + 8ci + 16c2 + 8c3 = 1. 



(2.4) 



We have to fix the gauge in perturbation calculation. We take the lattice Lorentz gauge given 
by 



(2.5) 



n fi 



Let us introduce the following Fourier transformation 



(2.6) 



where denotes n^=i Here and hereafter we make a change of variable from k^a to 



The free part of the action So is given by 
1 



where 



and 



So 



Tr q^u{k)F^u{k)F^A-k) + 2k^KA^{k)A,{-k) (2.7) 



k,, = 2 sin — . 

M 2 ' 



F^u — M kfiA/ylk) ki^A^(^k) j , 



(2.i 



(2.9) 



qf,u{k) = 1 - (ci - C2 - C3) yk^ + k^) - (c2 + C3) k ; for /i / z^, 
q^l^lik) = 0; for all fi. 



(2.10) 



Here 

4 

The free propagator is defined by 

>«(fe)yl^(fc')) = <5a6(2vr)4<5(4)(A; + A;')I?;..(fc), 
and is the inverse of the free two-point function 

G^yijv) = k^ky — ^ ^ {kp5^i, — k^5pi, \ q^pkp. 



(2.11) 



(2.12) 



(2.13) 



The explicit form of the propagator is given in the Appendix [Xj Note that the free part of 
action depends on only the combinations of ci — (c2 + C3) and C2 + C3. Therefore we have two 
parameters ci and C23 = 02 + 03 to choose in the lowest order perturbation theory, because cq is 
determined from them by the constraint (12. 4p . 



We will make the perturbation calculation with this free action. Let us denote the expectation 
value in the perturbation theory as 

(O) = J DAOe^'' I J DAe^°. (2.14) 

The free propagator of F^y{k) is defined by 

'F;u{k)Fl^{k')) = 5a,{2^f5^'\k + k')D,,,,^{k). (2.15) 

For example, D^y^p\{k) is given by 

D^,u,px{k) = (k^^^ D,, + (k,^^ D^^ - 2k^KD^,. (2.16) 

The expectation value of various Wilson loops can be written in the lowest order of pertur- 
bation theory as 

Ar2 - 1 

W{C) = l-g^^^F{C) (2.17) 

where C denotes the Wilson loop. Note that F{C) is independent of A^. The F{C) for some 
particular Wilson loops are given as follows: 

f /sin i/fci sin ijA;2\ 

F{I X J) = / 1)12,12 ■ \' . \. (2.18a) 
7a; y sm sm 2K2 j 

F(chair) = ^ L>i2,i2 Q + ^ cos k^ (2.18b) 
F(3-dim.) = ^ Z)i2,i2 Q + ^ cos k^ (2.18c) 
F(4-dim.) = / Dx2,\2 [6-2(1- cos k'i) (2 + cos k^)] (2.18d) 

Jk 

Here F{1 x J) is for / x J rectangular loop, i<'(chair) is for the loop depicted in Fig{T]c, F(3-dim.) 
for the loop in Fig|T]d. The 4-dimensional loop in eq. (12.18d|) may be described by the loop: 



n 



11+ i n + i + 2 n + i + 2 + 3 n + i + 2 + 3 + 4 
n + 2 + 3 + 4 n + 3 + 4 n + 4 n. 
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3. BLOCK SPIN TRANSFORMATION 



First note a difference between 2d 0(N) sigma model and 4d SU(N) gauge model in that the 
variables are site-variables for the former, while they are link-variables for the later. Therefore 
we propose to make blocks of link variables in 4d SU(N) gauge models, rather than blocks of 
sites. See FigH 

Secondly we are interested in the limit g — )• gc{= 0) and in the physical quantities for which 
perturbation theory is applicable. As discussed in the first section of [I], the correlation length 
is large near the critical point gc and diverges at the critical point as a result of cooperative 
behavior of the system. Within the correlation length the properties of the system do not change 
qualitatively. Therefore, if we make the weak coupling expansion, the variables A^{n) do not 
change rapidly within the correlation length. 



Thus we propose to define a block transformation by 
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Af,{2n) + YAf,{2n + i)) + YA^{2n + + p) + ^ A^{2n + i) + p + X) 

(3.1 



and 

C/(^) (n) = exp (igaA^^^ (n)) . (3.2) 



Note that the normalization constant is 1/8 rather than 1/16, because the length of the link 
in the new system is 2a. Except for the normalization, eq. ()3.ip is similar to eq.(3.1) in (I). 
Therefore let us write eq.(3.1) simply as 

4\n) = lY.Mn), (3.1a) 

ngn' 

and let us call this transformation block spin transformation, although the variable is not spin- 
variable. The block spin transformation is applied repeatedly. A block variable A^fj\n) after 
I — th iteration is defined as 

4'^(^') = ^E4'"'^H (3-3) 

nGn' 

and 

[/(^) (n) = exp (igaA^^^ (n)) . (3.4) 



FIG. 2. Schematic diagram for block spin transformation (3-diniensional projection). 



The effective action after the block spin transformation is defined through the relation 



e 



^«(^«) = / K f^(^-i), e'^'-''^^"-') (3.5) 



where 



M n' V nen' / 



(3.6) 



The renormalization group r is defined as 

r (5(^-1)) =S«. (3.7) 

The is the original action and y4^^^(n) is the original variable A^{n). This transformation 
is slightly different from that proposed by Wilson J], although it is the same in spirit. This 
difference becomes, however, crucial when ones consider the limit / 00, as will be shown in 
the next section. 

Let us calculate the expectation value of various Wilson loops for the I — th block link 
variables in the perturbation theory. Let us denote them as 

/V2 _ 1 

W{I){C) = l-g^^^F^'){C) (3.8) 

following eq. (l2.17p . The expectation value is taken with the action 5^^^ (yl^^-*). As is in 2d 0(N) 
sigma model, it is easy to calculate F^^^ (C) in terms of the original variables in the perturbation 
theory. We obtain 

F<')(/ X J) = / ( !i^£f (3.9a) 

\ sin 2^;) sin 2^ / 



Here 



f(^) (chair) = I Dg\, Q + ^cosA;f ) H^'\k), (3.9b) 
F(^)(3-dim) = I D% + ^cosfcf ) H^'\k), (3.9c) 
F(^)(4-dim) = J dI^Iu 6- 2(^1 -cos 4^^) (^2 + cos 4^^) H^^\k). (3.9d) 



k^P = 2'k^ (/i = l~4) (3.10) 

/-I 4 



4 

M=0 Ai=l 



and 



^iiU(^) = {k\i^)^ D,,{k) + (4^))'Z)^^(A;) 

- 2fe(.^)4^) cos (^(^2^-1 - k^ cos (^(^2^-1 - k,^ D^,{k) (3.12) 



where 



4^)=2sin^. (3.13) 

The derivation of eq. (j3.9p is given in the Appendix |Bl 

The gauge for the block variables (n) {I ^ 1) is not the lattice Lorentz gauge, in general. 
This is not the problem to calculate F^^\C), as is shown above. 



4. EXPECTATION VALUE OF VARIOUS WILSON LOOPS ON THE 
RENORMALIZED TRAJECTORY 

In the previous section we have obtained the expectation value of various Wilson loops for 
block link variables after I — th iteration. In the limit / 00, the effective action S^^^ should 
approach the renormalized trajectory. Therefore the expectation value of the Wilson loop also 
should approach to that on the renormalized trajectory. Let us denote 

lim W'^^XC) = W^^\C) = 1 - g^-^-T7^F^^\C) (4.1) 
/— )-oo 4A' 



In this section we will derive 



First note 



lim / e'^^'D^.ik) 

n— >-oo J 



47r2 n2 '^'^ 



(4.2) 



The first term in the r.h.s. is the dominant contribution to the expectation value in the limit 
/ — )• oo. As in the section 4 of [I], we are able to prove that the non-leading terms which come 
from the second term of the r.h.s. of eq. (l4.2p do not contribute to the expectation value in the 
limit I — )• oo. Hence we neglect them. Therefore we have to consider only D^^{k) with fi = ly. 
Secondly let us remind that 



2^-1 



,m,n,T) = 



(4.3) 



Zl ,mi 
ni,ri 



Now let us consider, for example, ^^^-'(1 x 1). From eqs. ()4.2p and (j4.3p we have 
f(^)(1 X 1) ~ 4 f ^ j [F{li-l2,mi-m2,ni-n2,Ti-T2) 

-F[2^ + h - h^rni — m2,rai - ra2,Ti - T2) 



h,m\,l2,m2 
ni,T\,n2,T2 



where 



F(/, m, n, r) 



1 



(4.4) 



(4.5) 



The approximate equality in eq. (|4.4p means that only the leading terms are taken. Equation 
(|4.4p may be written as 



1 



Tl — T2 ' 

21 ; 



1 + ^) +(^^v^)' + (^v^)' + (^)'. 



(4.6) 



Note that the factor (1/2)^^ are multiplied for the denominator and the coefficient. Putting 
1/2^ = e as in [I], we obtain finally 



^(00) X 1) = y dx\dy\dz\dT\ J dx2dy2dz2dT2 



(Xi - X2f + (yi - ^2)^ + {Zl - Z2f + (ti - T2)^ 



1 



(1 + xi - X2) + (yi - 7/2) + {zi - Z2) + (n - T2) 



(4.7) 



9 



Introducing the function F defined by 

1 



F(ni,n2,n3,n4) 



E 



1 



47r2 ,j^^f;^3^,4 (ni - eix)^ + (n2 - £2^)^ + (ns - esz)^ + (n4 - e^rf 
±1 



(4.8) 



and setting xi — X2 = x and so on, we obtain 

1 

F(°°)(1x1)=4 jjjj (ixdydzdr(l-x)(l-y)(l-z)(l-r)x (|f(0, 0,0,0) 0,0,0)) . (4.9) 


We are able to obtain other F^°°\C) similarly. 



Let us introduce the symbol 



(1 — x){l — y){l — z){l — T)dxdydzdT 

for notational simplicity. Then we have, for example, 

fM(1 xl) = 4[jj] (f(0, 0,0,0) -#(1,0,0,0)) 

F^°°\l X 2) = 6F(0, 0,0,0) -4F(1, 1,0,0) -F(2, 0,0,0) 

f(°°)(2 X 2) = F(0,0,0,0) +F(1, 0,0,0) -F(2, 0,0,0) -F(2, 1,0,0) 



(4.10) 



(chair) = 6F(0, 0,0,0) -4F(1, 0,0,0) -2F(1, 1,0,0) 
F(°°)(3-dim) = 6 F(0, 0,0,0) -F(l, 1,0,0) 



F(°°)(4-dim) 



F(0,0,0,0) -F(l, 1,1,0) 



(4.11) 



It should be noted that it is non-trivial that F^^^C) approach finite constant value in the 
limit / —7- 00. It depends on the definition of renormalization group. If we had taken the 
renormalization group proposed in the ref.[4], we would obtain that F^°°\C) are trivially zero: 
Because the length of the new link after the renormalization does not change from that of the 
original link, we have the sum over only li, mi, and ni in the equation corresponding to eq. ()4.3p 
and consequently we have an extra factor (1/2^)^ which vanishes in the limit / 00, in the 
equation corresponding to eq. (l4.7p . 



5. BLOCK SPIN TRANSFORMATION FOR VARIOUS ACTIONS 



Now we are ready to calculate the function F^^^ (C) of the block variables for any parameter 
ci and C23 (= C2 + C3), using eqs. (|3.9|) ~ (|3.13p . We first list some of them for the standard model 
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(ci = C2 = C3 = 0) together with F^°°\C) obtained from eqs. (j4.10p and (|4.1ip . We clearly see 
from the table |T] that the functions gradually approach the asymptotic values. We also list some 
of F^^\C) in the table |T] for ci = — 0.252, C2 = 0, C3 = —0.17 (model W) which has been chosen 

yy Wilson in ref.|^, for ci = — 1/12,C2 = 03 = (model WZ) which has been chosen by Weisz 

I 1 

3| in Symanzik approach [5], as well as for ci = — 0.331, C2 = C3 = (model IMll) and for 
ci = —0.27, C23 = —0.04 (model IM22) which will be chosen by a certain criterion below. In 
Figj3]we depict block spin transformation flows for the function F^^\C). 

We see from the Figure and the Table that the behavior of the convergence for the function 
to the fixed point crucially on the parameter ci and C23. 



6. IMPROVED LATTICE ACTION 

According to our strategy described in detail in [I], let us choose an action which is located 
near the renormalized trajectory. First let us note that there is no difference between interac- 
tion, for example, Tr {U1U2U3U4) Tr {U^^U^UqUj) and Tr {U1U2U3U5UQU7) in the lowest order 
perturbation theory. Therefore we assume that the action is of the form 

S(°^) = ^KcW{C) (6.1) 
c 

where W{C) is the Wilson loop for the contour C. Practically we have to truncate the sum in 
eq. (l6.ip . Let us restrict actions to those given by eq. (|2.ip . 

We define a distance from an action to the renormalized trajectory by 



where M is the number of terms in the sum over the contour C. We restrict the contour in the 
sum to those up to six-link length. 

Plotting and R^'^^ defined above in the two-dimensional space spanned by ci and C23, we 
find that there is a one-dimensional very narrow deep valley where ii^^^ and R^'^^ are very small, 
respectively. For R^^^ , the one-dimensional line is parameterized by the equation 

1.531C1 + C23 = -0.5067 (6.3) 

and for 

1.746C1 + C23 = -0.5116. (6.4) 
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TABLE I. The values ol F(^)(C) for various models together with i^(°°)(C). 
Table 1-a Model S (ci = ca = C3 = 0) 





r (^i X i j 


r [1 X Z) 


r [Z X Z) 


r (^cnaii J 


r yo-dim.) 


r (^4-aim.j 




0.50000 


0.86225 


1.36931 


0.78444 


0.85331 


1.17759 




0.28810 


0.51765 


0.87978 


0.45667 


0.50569 


70854 


F(2) 


0.21623 


0.40352 


0.72087 


0.34655 


0.39094 


0.55476 


i^(3) 


0.19446 


0.36985 


0.67515 


0.31350 


0.35711 


0.50956 


i^(4) 


0.18865 


0.36183 


0.66877 


0.30465 


0.34801 


0.49680 




0.18649 


0.35770 


0.65875 


0.30146 


0.34493 


49331 


Table 1-b 


Model W (ci = 


-0.252, C2 0,C3 = 


-0.17) 




F(l X 1) 


F(l X 2) 


F(2 X 2) 


F (chair) 


F(3-dim.) 


F(4-dim.) 




0.19297 


0.36601 


0.64859 


0.31059 


0.35287 


50047 




0.18211 


0.34432 


0.62927 


0.29439 


0.33685 


0.48168 


i^(2) 


0.18271 


0.34871 


0.64290 


0.29573 


0.33906 


0.48554 




0.18649 


0.35770 


0.65875 


0.30146 


0.34493 


0.49331 




Table 1-c 


Model WZ (ci = -1/12, C2 = cg 


— n^i 

— ^) 




F(l X 1) 


F(l X 2) 


F(2 X 2) 


F (chair) 


F(3-dim.) 


r (^4-aim.J 


F(o) 


0.36626 


0.66263 


1.09814 


0.57705 


0.63235 


0.87786 




0.25076 


0.46043 


0.80091 


0.39900 


0.44470 


0.62629 


F(2) 


0.20599 


0.38751 


0.69875 


0.33087 


0.37464 


0.53292 


F(3) 


0.19173 


0.36556 


0.66924 


0.30934 


0.35285 


0.50387 


F(4) 


0.18794 


0.36072 


0.66724 


0.30358 


0.34691 


0.49534 




0.18649 


0.35770 


0.65875 


0.30146 


0.34493 


0.49331 


Table 1-d 


Model IMll ( 


ci ^ -0.331, C2 = C3 = 0) 




F(l X 1) 


F(l X 2) 


F(2 X 2) 


F (chair) 


F(3-dim.) 


F(4-dim.) 


i?(o) 


0.21027 


0.40340 


0.71109 


0.33352 


0.36977 


0.51874 




0.18826 


0.35756 


0.65039 


0.30184 


0.34074 


0.48500 


F(2) 


0.18431 


0.35247 


0.64945 


0.29743 


0.33935 


0.48528 


j!?(oo) 


0.18649 


0.35770 


0.65875 


0.30146 


0.34493 


0.49331 



Along the line defined by eq. (|6.3p . R^^^ is less than 0.01 and the variation of R^^'^ is not so 
rapid as far as > C23 ^ —0.2. The minimum of i?*-^) is about 0.00545 at ci = —0.27 and 
C23 = —0.09333 (to be referred as point IM12). When we put C23 = 0, the minimum of R^^^ is 
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TABLE I. (Continue) 
Table 1-e Model IM22 (ci ^ -0.27, 023 = -0.04) 





F(l X 1) 


F(l X 2) 


F(2 X 2) 


F (chair) 


F(3-dim.) 


F(4-dim.) 


F(0) 
i^(2) 


0.22081 
0.19423 
0.18701 


0.41968 
0.36707 
0.35671 


0.73586 
0.66472 
0.65556 


0.35136 
0.31177 
0.30181 


0.39167 
0.35263 
0.34439 


0.55003 
0.50167 
0.49220 




0.18649 


0.35770 


0.65875 


0.30146 


0.34493 


0.49331 



about 0.00773 at ci = —0.331 (point IMll). Note that the values of R^^^ do not differ so much 
between the two points. 

Along the line defined by eq. (j6.4|) . is less than 0.0025 and the variation of R^'^^ is not 
so rapid as far as > C23 > —0.2. The minimum of i?^^^ is about 0.00219 at ci = —0.27 and 

= -0.04 (point IM22). When we put C23 = 0, the minimum is about 0.00242 at ci = -0.293 
(point IM21). 

If we could perform MC simulations with parameters ci and C2 (or C3) on a very large lattice, 
the action IM22 could be the best action among those considered up to here. However, it makes 
a large difference for computer time in MC simulations whether we include the C2 (or C3) term 
or not. On the other hand there is no large difference between action IM22 and action IM21, or 
between action IM12 and action IMll from the view point of renormalization group. Therefore 
it is practically better to choose action IM21 or action IMll. 

There is also a limitation on the size of the lattice where we make MC simulations. For 
example, we will measure the string tension 6( on a 8^ lattice for SU(3) lattice gauge theory. In 
this case the string tension is mainly determined by the Creutz ratio [7] X{3,3) (and X(4, 4)). 
The X{3, 3) is determined from W{3 x 3), VF(2 x 3) and W{2 x 2). On the other hand, ^(^^(l x 2), 
for example, contains the contributions from the propagator Di^i^{m,n) up to m = 3 and n = 5, 
while ^'^^^(l X 2) contains those up to m = 7 and n = 11. Thus to obtain the Wilson loops 
such as W{3 x 3) or W{4 x 4), the criterion based on R^^^ is better than on R^'^\ From this 
consideration and an analysis of instantons on the lattice (see section [8]) , we choose the model 
IMll as an improved action (model IM). When we make MC simulations on a lattice, e. g., 16^, 
we may choose the model IM21 as an improved action. Anyway the difference concerning R(n 
between models IMll and IM21 is not large. 



In ref.ja], we have measured the string tension on a 8 lattice for the SU(2) lattice gauge 

13 



u 


. y 












AO 





.8 

















. 7 








■ 









.6 

















.5 


- 




Ai 

■ 1 


A 




- 





. 4 






, — ^ Fig . 3-b 


model 


s 












■ 

T 


model 


wz 





. 3 








model 


w 


1 

Fixed Point 




• 


model 


IM22 





.2 



















. 1 


0.2 


0.3 


0.4 




0.5 



F (IXl) 



0.4 



0.38 



CM 

X 



0.36 



0.34 



0.32 




0.17 



. 18 



.21 



.22 



(b) 

FIG. 3. Block spin renormalization group flows for X 1) and F'-^\1 x 2) for various models. The 
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14 



theory with the action ci = —0.3371, C2 = C3 = (model R3). This action has been chosen 
from the consideration of the scale parameter and an analysis of instantons on the lattice. For 
this action R^^'^ = 0.0097. Compare this value with i?^^) = 0.027 for model W, = 0.312 
for model WZ and R^^^ = 0.495 for model S. Thus the model R3 is also an improved action 
compared with models WZ and S: As far as ci ~ —0.3 (with C23 = 0) for which R^^^ < 0.01, any 
action can be taken as an improved action. 

In ref. 8| we have chosen In Al/A^ as a parameter for the distance between an action and 
action W. It seems that a line where R^'^\ for example, is constant corresponds approximately 
to a line where the scale parameter A is constant, as noted in [I] for 2d 0(N) sigma models. 
Thus the criterion chosen in ref.[^] for the distance is reasonable even from the view point of 
the renormalization group. 

Note that model WZ is not an improved action due to our criteria: The value R^^^ = 0.312 
implies that model WZ is far from the renormalization trajectory, and the life time of instantons 
on the lattice is short. On the other hand, model W is an improved action; although R^^^ is 
slightly larger than that of model IMll or that of model R3, it is of the same order. 



7. IMPLICATIONS FOR MC CALCULATIONS 



As far we have already discussed the general feature of the renormalization group and our 
strategy in section 7 of [I], we do not repeat them here. According to our strategy, let us set an 
upper limit A for the relative difference between F^^\C) and F^^\C) as A = 0.01. Then for 
model IM we have J = 1 as the minimum number of iteration for which i?*^"^) < A is satisfied. 
On the other hand for model S we have J = 4. 

This implies that the Wilson loops W{I x J) agree with W^°°\l x J) within relative difference 
0.01 for /, J > 2 in the case of model IM, according to our discussion in section 7 of [I]. This 
further implies that the Creutz ratio A(/, J) take their asymptotic value on the renormalized 
trajectory for I,J>3. 

On the other hand, the Wilson loop W{I x J) for model IM corresponds to the Wilson loop 
W (2^1 X 2^ j) for model S. Thus it is only expected that the Creutz ratios X{I, J) take their 
asymptotic values on the renormalization trajectory, for /, J > 16 with precision 0.01 in the case 
of model S. (See FigH]) This is the reason why we call IM the improved action. 
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FIG. 4. Schematic diagram of the renormahzation group from model S to model IM. The model IM in the 
Figure does not exactly correspond to the model IM defined in the text. However, they are approximately 
equivalent to each other, because i?^^) for both models are of the same other. 



8. INSTANTONS ON THE LATTICE AND THE RENORMALIZED TRAJECTORY 

The general discussion on instantons on the lattice given in section 10 of [I] can be also 
applied for 4d SU(N) gauge models. Thus we conclude also for 4d SU(N) gauge models that the 
renormalized trajectory is located at the boundary which devides the parameter space into two 
parts: In one of them instantons exist, while in the other instantons do not exist. We further 
conclude that the one-dimensional line such as defined by eq. ()6.4p divides the two-dimensional 
space spanned by ci and C23 into the two parts. See FigiH 

It is rather difficult to verify this conclusion numerically compared with the case of 2d 0(3) 
sigma model, because it takes a lot of computer time to do it. Therefore we have investigated 
y, 0] the existence of instantons on a lattice with size of 6^ with ci being varied and with 
C2 = C3 = by the method described in ref.jol. We have found the following: For model 
R3 (ci = —0.3375) we have instantons with topological number q = 1,2,3,4 for ten random 
starts, for model IMll (ci = —0.331) we have instantons with q = 1,2,3 and for model IM21 
(ci = —0.293) we have instantons with q = 2,3,4. For |ci| < 0.29, we have no stable instantons 
on the lattice. Thus the point ci = —0.29 is critical. This is consistent with the above conclusion. 

One reason in addition to the reason given in section [6] why we choose model IMll rather 
than model IM21 as an improved action for a 8^ lattice is that we have no instantons with q = 1 
for model IM21 on a 6^ lattice as far as we have investigated. We expect that if the size of the 
lattice is large enough, e.g., 16^, instantons with q = 1 also will exist on the lattice for model 
IM21. 
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FIG. 5. Stability of instantons on the lattice vs. ci and C23. 



9. DISCUSSION 



The general discussion given in section 11 of [I] can be also applied for 4d SU(N) gauge 
models. We state here only briefly the main points. The improved action near the renormalized 
trajectory has been determined by the perturbation theory in our approach. We are able to 
calculate non-perturbative effects using this action by MC simulations. The effect of instantons, 
for example, is properly taken into account even on a lattice with small size. 

A different approach to improve lattice action is proposed by Symanzik Q]. The perturbation 
theory cannot be applied for calculation of physical quantities such as the string tension where 
non-perturbative effects are crucial. It appears that this limitation is not properly taken into 
consideration in Symazik's approach. As shown in section [6l the "improved" action obtained by 
Weisz [3] in Symanzik's approach is not an improved action in our criteria. 

The continuum limit of a lattice theory is unique and universal for a wide class of lattice 
action because of the following facts: i) The renormalized trajectory is unique and ii) the 
trajectory for any lattice action approaches asymptotically the renormalized trajectory, if the 
action belongs to the domain in the parameter space which is governed by the non-trivial fixed 
point. However, it crucially depends on the form of the action how rapidly the trajectory 
approaches the renormalized trajectory. We have given a criterion to estimate this rapidity. 
From the criterion we can guess the coupling constant where we may expect that the expectation 
value of a physical quantity is identical with that on the renormalized trajectory. The point is 
not that a physical quantity approximately shows a scaling form, but that it is identical with 
that on the renormalized trajectory. 

We would like to calculate various physical quantities with an improved action due to our 
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criteria in addition to the string tension of the SU(3) gauge theory [g]. We hope that the mass 
ratios of mesons to baryons will become realistic with an improved action of a relatively small 
lattice. We also conjecture that the first order phase transition [10] observed in SU(4) and SU(5) 
gauge theories with the standard action will not be observed with an improved action, because 
no phase transition will occur on the renormalized trajectory. 



We finally argue that Osterwalder-Schrader 



111 ] (OS) positivity is satisfied on the renormal- 



ized trajectory and consequently it is satisfied, at least approximately, for the improved action 
of 4d SU(N) gauge models, from the same reasoning given in [I]. 
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Appendix A: The explicit form of the propagator [3| 



It is easy to show that Dfj_^{k) may be written in the form 



-2 



k^kp + ^ ^ {kcj5^i, kpd^uj -^uaka 



(A.l) 



with Af^i, satisfying (i) = for ah fi and (ii) A^,^ = A^^. 

The element A12 (and the other elements by appropriate replacement of indices) is given by 



A 



12 



- ^2 ) (913914^? + qisqsAkl + qms^kj 



+ [k -k^] ( 923924^:2 + 923934^:3 + Q'24g34A:4 





+ gi3g24 (kf + fci) (k^ + kjj + ^14923 (kf + kjj (kl + kl 



^12^34 [k^ + k^j - (^13923 + 914924) ^3^4 

912 [qiskjkj + qukjki + 923^2^4 + 924^1^1 



(A.2) 



where 



A 



XI n ^^'^ + X] k^klqi^y {q^^pqur + q,iTqup)- 



(A.3) 



Appendix B: Derivation of ea. (|3.9p 



Let us introduce the following Fourier transformation 



4^)(n)= / e^2^('="+^'^/2)4^)(A:), 



(B.l) 

Jk 

because the system for variables A^fj\n) is that with spacing 2^a. Combining eqs. (|2.6p . ()3.3p 
and (jB.ip we obtain 



4)(A;) = e-'i''""''/^>-m'\k)A'-^\k), 



(B.2) 



where 



7-1 , 4 



(B.3) 



M=0 M=l 



19 



Let us define the free propagator L>(^) for the field A^^) by 

{Af\k)Al^^\k')) = 5ab{27r)^5(^\k + k^jD^jJik). (B.4) 
Then from eq. ()B.2p we have 

D(!J{k) = e-<''~'~'/'>-e~<'''"-'/^>^H('Hk)D<i^]{k), (B.5) 
where H'-^^k) is defined by ea. (l3TT]l and D'jj^Jik) = D^,^{k). 

The expectation value of various Wilson loops may be written as 

W{C) = J2(^f^Ak)D^,Ak) (B.6) 

for the original system. Then the corresponding expectation value for the I — th system can be 
written as 

W^'\C) = Y,c,A2'k)Di!J{k), (B.7) 

because the I — th system is that which scaling 2^a. Note that D^u{k) is odd in k^ and ky when 
^ / z^. Thus from eqs. flB.Sp and (|B.7p we can derive eq. (l3.9p . 
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